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Abstract 

Entanglement distillation is a procedure for extracting from one or more pairs of entangled 
qubits a smaller number of pairs with a higher degree of entanglement that is essential for many 
applications in quantum information science. Schmidt projection is a unique distillation method 
that can be applied to more than two pairs of entangled qubits at the same time with high efficiency. 
We implement the Schmidt projection protocol by applying single-photon two-qubit quantum logic 
to a partially hyperentangled four-qubit input state. The basic elements of the procedure are 
characterized and we confirm that the output of Schmidt projection is always maximally entangled 
independent of the input entanglement quality and that the measured distillation efficiency agrees 
with theoretical predictions. We show that our implementation strategy can be applied to other 
quantum systems such as trapped atoms and trapped ions. 
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I. INTRODUCTION 



Shared entanglement is an essential resource in many quantum information processing 
applications such as quantum teleportation quantum computation , quantum 

cryptography j^, and quantum repeaters j^. These applications generally require that two 
remote parties maintain a very high degree of entanglement in order for their protocols to 
work properly. Therefore, the generation of maximally entangled states such as the Bell 
states is an essential task in quantum information science. However, even if one can initially 
generate maximally entangled states with ease, they often become non-maximally entangled 
or partially mixed at the remote locations because of dissipation and decoherence through 
interactions with the environment. While each application has its own level of tolerance to 
entanglement degradation, most applications would be practically unusable without some 
kind of procedure to restore the degraded entanglement to its initial maximally entangled 
state. 

The need for most quantum information processing applications to overcome unavoidable 
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entanglement has led to several methods for entanglement restoration 



0, y, y, Q, 



13l . 1 14 Il5l |. These techniques typically seek to improve the quality of entanglement 
between a pair of entangled qubits by selecting, filtering, and applying unique processing 
algorithms to a larger number of lower quality entangled qubits. Entanglement restoration 
provides a path to utilizing entangled qubits that would otherwise be dicarded because 
their entanglement quality is below the threshold necessary for specific applications. The 
cost of entanglement restoration is that many pairs of the lower quality entangled qubits 
must be generated and processed, which requires more resources and longer times for each 
application processing step. It is therefore of interest to choose a method that is efficient 
in terms of resources and processing required to achieve a certain level of entanglement 
quality. Among these methods, Schmidt projection (SP) Q] has the interesting property 
that it can be applied to more than two pairs of entangled qubits at the same time and that 
the distillation efficiency can approach unity when applied to a large number of initial pairs. 
Despite these advantages, Schmidt projection is considered difficult to implement because 
it requires simultaneous collective measurements on multiple qubits. ^ 

In this work, we illustrate the working principle of Schmidt projection [7] by applying the 
technique to two pairs of partially entangled qubits. First, we use the concept of physical 
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mapping to encode the two pairs of entangled qubits into four physically distinguishable 
states of two hyperentangled photons. The two photons are entangled in both polarization 
and momentum degrees of freedom so that the four qubits can be distinguished by their pho- 
ton, polarization, and momentum. By hosting two qubits per photon carrier, we eliminate 
the problem of bringing the two qubits together to overlap completely in time for quantum 
logic gate operations. Moreover, the two types of qubits can be manipulated determinis- 
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tically using single-photon two-qubit (SPTQ) quantum logic 
carried out efficiently using SPTQ quantum gates to produce a single maximally-entangled 
pair of qubits. In Sect. II we describe the basic idea of the SP protocol and show how this 
protocol can be implemented with hyperentangled photon pairs. Section III describes our 
experimental implementation including the generation of partially hyperentangled photons 
and the use of SPTQ logic to perform the Schmidt projection. The resulting maximally 
polarization-entangled photons are characterized and the effectiveness and efficiency of the 
protocol are analyzed and discussed in Sect. IV, and finally we summarize and show how 
the concept of physical mapping can be applied to trapped-atom and trapped-ion systems. 

II. THEORETICAL BACKGROUND 

Before we describe how the Schmidt projection protocol can be implemented in our experi- 
ment, it is useful to clarify the sometimes confusing terminology in entanglement restoration. 
Given one or more pairs of entangled qubits one can extract from them a smaller number 
of pairs with a higher degree of entanglement (or the same number of pairs with less than 
unity probability) by using local operations and classical communication. This process of 
improving the amount of entanglement is called entanglement distillation. In the literature, 
this process is sometimes called entanglement purification or entanglement concentration. 



In this work, we follow the convention given by Ref. [lOl], in which entanglement purification 
refers to the process of enhancing the purity of a mixed state by extracting a number of 
purer entangled states from a larger number of less pure entangled pairs. Entanglement 
concentration is the process for increasing both the degree of entanglement and the purity 
of a given initial state. 

To be more specific, we use the definition suggested by Bennett et al. [7,] to quantify 
entanglement. For a partially entangled pure state shared by Alice {A) and Bob (5), 
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the von Neumann entropy of the partial density matrix gives a measure of the entanglement 

E{ip) = -Tr(p^log2PA) = -Tr(pijlog2Pij) , (1) 
where Pa{pb) is the partial trace of \ip){il)\ over subsystem B {A). In general, the en- 
tanglement measure is bounded by < E{;ip) < log2 d for a bipartite qudit system in a 
2(i-dimensional Hilbert space. A natural extension of the definition to a mixed state with a 
density matrix p is the entanglement of formation E{p), which is the smallest expectation 
value of entanglement E of any ensemble of pure states realizing p IGj. In other words, to 
calculate the entanglement of formation, one needs to consider all possible decompositions of 
p, that is, all ensembles of states with probabilities Pi satisfying p = Sjpj|^/'j) Then 
the entanglement of formation is defined as E{p) = mmTiiPiE{tpi), which in general is diffi- 
cult to calculate. However, for a bipartite qubit system, an explicit formula for entanglement 
of formation exists llj 



Most of the protocols for enhancin g E can only be applied to one |10l . |12| ] or two pairs 
of entangled qubits at a time 14. Iisl and therefore have a fixed distillation efficiency 
independent of the number of initial pairs. In contrast, Bennett et al. showed that Schmidt 
projection can be applied to n pairs of partially entangled pure state with E{;ip) to extract 
n ■ E{%1)) pairs of maximally entangled Bell state in the limit of large n [3], and this property 
is used to justify the definition of E in Eq. ([1]) as a measure of the degree of entanglement. 
Entanglement Eltp), and hence the SP technique, forms the basis of or is related to other 



common measures such as entanglement of formation, concurrence [18|, ll9|, and tangle |20 |. 
Despite its role in the theoretical foundation of quantum information science, Schmidt pro- 
jection has not been experimentally demonstrated due to the difficulty of performing the 
required collective measurements. 

To understand how Schmidt projection works, consider two pairs of identically entangled 
state, \4')ab = cos6'|0)a ®\0)b + sin 6*11)^ ® which can be described by 

= cos^ ^|0)^JO)aJO)bJO)b, + cos e sin ^|0)aJ1)a2 |0)bJ1)b2 

+ sinecose|l)AjO)^,|l)BjO)B2 + sin2^|l)Ajl)Ajl)i?Jl)B2, (2) 

= cos^9\0)a\0)b + sin2^|3)^|3)B + cos ^ sin ^(| 1)^1 + \2)a\2)b) • (3) 

In Eq. ([3]) we replace the binary representation of Alice's qubits and Bob's qubits by qudit 



notation (decimal number). For two entangled pairs, the task of collective measurements is 
reduced to projecting terms with the same coefficient (cos ^ sin whose output form is that 
of a maximally entangled state. Equation shows that each of the three terms is locally 
orthogonal to each other, and therefore the result of Alice's local projection is perfectly 
correlated to that of Bob's local projection, thus making post-processing and the associated 
classical communication unnecessary. In contrast, the non-Schmidt projection methods used 
in Ref . require that Alice and Bob compare their measurement results over a classical 

channel and change the phase based on their measurements. 

Schmidt projection is most useful when it is applied to n > 2 pairs of \4')ab that has the 
initial tensor product state given by 

mi 

= (cos^|0)aJO)b, + sme\l)AAl)B,) ® ■ ■ ■ ® (cos0|O)aJO)b„ + sin^|l)^Jl)Bj 
= cos"^|OU...|0)AjO)B,...|0)i,„ 

+ cos"-^0sin0(|O)A, . . . |0)^„_J1)aJ0)b, • • • |0)b„_J1)b„ + ■ ■ • 

+ \1)aA0)a,...\0)aJ1)bA0)b,...\0)bJ 
-f--- + sin"e|l)^,...|l)^Jl)B,...|lV. (4) 

Among the 2" terms in Eq. (jlj), there are n + 1 distinct Schmidt coefficients (cos"^, 
cos"~^ 9 sin 9, sin" 9), and each group with the same coefficient cos"~'^ 9 sin^ 9 has (^) terms 
that are maximally entangled, for n > k > [7]. Alice and Bob project this initial state 
into one of — 1 orthogonal subspaces to extract the associated maximally-entangled states. 
The extracted state can be directly used for faithful teleportation in a (^)-dimensional or 
smaller Hilbert space, or, alternatively, one can efficiently convert this state into a tensor 
product of Bell states, as suggested by Bennett et al. 

In our experiment, we demonstrate a proof-of-principle implementation of Schmidt pro- 
jection by applying the SP protocol to two pairs of entangled qubits such as those given 
by Eq. ([3]). We choose to implement the four-qubit SP protocol with a pair of partially 
hyperentangled photons as the initial state given by 

l^p) ® I^m) = {cos9p\V)a\V)b + ^in9p\H)A\H)B) 

®{cos9m\L)a\L)b + ^i^9m\R)a\R)b) , (5) 
where subscript P (M) refers to the polarization (momentum) degree of freedom, H iV) 



represents horizontal (vertical) polarization, and L (R) is the left (right) path shown in 
Fig. [TJ Schmidt projection requires that the two input qubit pairs are identically entangled, 
Op = 9m = 6, which we are able to set by the method described in the next section. With 
this setting, Eq. can be expanded as 

cos^ e\VL)A\VL)B + sin^ e\HR)A\HR)B 

+ cose sin 9 {\VR)a\VR)b + \HL)a\HL)b). (6) 

In Eq. the SP protocol is applied by projecting the term associated with the coefficient 
cos 6 sin 6 that can be realized by transmitting the ^-polarized photon in the R path and 
the if-polarized photon in the L path by both Alice and Bob. 



III. EXPERIMENTAL SETUP 



The use of hyperentangled photons for demonstrating a four-qubit SP protocol has a 
number of experimental advantages. Only a pair of hyperentangled photons is needed to 
produce two pairs of entangled qubits, one pair entangled in the polarization mode and the 
other pair in the momentum (or path) mode. Therefore, the photon serves as a carrier 
for the polarization and momentum qubits that can be manipulated with single-qubit and 
;wo-qubit quantum gates using deterministic single- photon two-qubit quantum logic gates 
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22l |. We have previously demonstrated deterministic SPTQ quantum logic gates such as 



polarization- CNOT [21], momentum-CNOT and SWAP [22| using linear optical components. 
The SPTQ implementation is simple and offers a useful platform for few-qubit quantum in- 
formation processing tasks. In this section, we describe the method for generating partially 
hyperentangled photon pairs with adjustable degrees of entanglement using a bidirectionally 
pumped polarization Sagnac interferometer (PSI). We then detail the experimental imple- 
mentation of the SP protocol for extracting a pair of maximally-entangled qubits from the 
initial partially hyperentangled four-qubit state using SPTQ quantum logic. The setup for 
verifying the input and output states of the distillation procedure will be described. 



A. Generation of partially hyperentangled photon pairs 

To generate the hyperentangled state that are entangled in both polarization and momen- 
tum degrees of freedom, we utilized a polarization-entangled photon pair source based on 
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Idler (To Bob) 




FIG. 1: (Color online) Hyperentangled photon pair source in a PSI configuration, showing gener- 
ation of polarization-entangled photon pairs in (a) collinear output modes along Ra and Rb and 
(b) non-collinear output modes along La and Lb- The PSI consists of the crystal, two flat mirrors, 
a dual- wavelength half- wave plate (HWP), and a dual- wavelength polarizing beam splitter (PBS). 
The hyperentangled output state is a superposition of these two output modes each with a different 
probability amplitude. The two HWPs are used to rotate the polarizations by 90°. The subscripts 
1 and 2 for the signal and idler outputs in (b) refer to the clockwise and counter-clockwise directions 
of the SPDC outputs, respectively. DM: dichroic mirror, DA: double aperture mask. 

spontaneous parametric down-conversion (SPDC) in a type-II phase-matched periodi cally 



poled KTiOP04 (PPKTP) crystal in a bidirectionally pumped PSI configuration 23|, |24| . 
The SPDC source was driven with ~5 mW of a continuous- wave (cw) diode laser operating 
at a wavelength of 404.775 nm that was delivered to the PSI in a single-mode optical fiber. 
The 1-cm long PPKTP crystal was set at a temperature of ~19°C for nearly frequency- 
degenerate operation. The signal and idler outputs at 809.55±0.5 nm were detected through 
1-nm interference filters (IF) to restrict the SPDC output bandwidth. More details of the 



PSI source can be found in Ref. 
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\2M- 

Normally the PSI entanglement source is operated in a coUinearly propagating geometry 
for high collection efficiency of the polarization-entangled photons. However, it can also 
be operated in a slightly non-collinear configuration to exploit the momentum correlation 
of SPDC outputs. By collecting the non-collinear output modes from the PPKTP crystal 



the output photons are also momentum entangled 
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251], as illustrated schematically in 
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Fig. [H After the PSI outputs are collimated the momentum of the output mode can be iden- 
tified equivalently by its path location that is defined by the double-aperture (DA) mask. 



Unlike other hyperentanglement sources 
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28| , our experiment requires partial hyper- 



entanglement with adjustable but equal amounts of entanglement in both polarization and 
momentum degrees of freedom, as indicated in Eq. To control the degree of momentum 
entanglement (Om), we aligned the DA mask at an asymmetric location with respect to the 
pump beam axis so that the right aperture {Ra or Rb) coincided with the pump axis as 
shown in Fig. [H The diameter of each aperture of the DA mask was 1 mm, and the distance 
between the centers of two apertures was 2 mm. 

Consider the collinear R^ and Rb output modes of Fig. [I](a), which is the standard PSI 
configuration for generating polarization-entangled photons. We showed in previous PSI 



experiments 
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2J] that the output state is proportional to 

|V^)fiAi?fl = ((^osOpa^Hj^J'lj^^ + e"^ sine pal j^J'^jjjiJlO) , (7) 



where o,y j^^(b^jj j^^) is a (H) polarized photon creation operator for Ahce's Ra (Bob's Rb) 
path, and similarly for a|/^^ and b^n Ra- '^^^ relative amplitude {9p) and phase (0) of the 
output state are determined by the adjustable relative amplitude and phase between the H 
and V polarization components of the pump. Similarly, for the non-collinear La-Lb output 
modes of Fig. [I](b), the output state is proportional to 

\i^)LALB = (cos^pal^^^Jt.^^^ + sin ^pat.^^ Jl,^^ J 1 0) , (8) 

where 6p and (j) are the same as those in Eq. ([7]) because they are generated by the same 
pump laser. 

For hyperentanglement, we consider the coherently driven outputs of all four paths and 
obtain a final state that is a superposition of the \'ip)R^Rg and LaLb outputs of Eq. ([7]) and 
Eq. ([8]), respectively. In general, the collinear output modes (a^^, b^Rg) and the non-collinear 
output modes ^\b^ excited with different probability amplitudes depending on the 
phase- matched emission angle [29!]. Therefore, the four-path output is a partially hyper- 
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entangled state given by 



CO?, 6 m\'iP)laLb + sin 6'm I V')RARfl 
= cos9M{cos9p\HL)A\VL)B + e'^sinep\VL)A\HL)B) 
+ ?>m9M{cosep\HR)A\VR)B + e'^ ?>mep\VR)A\HR)B) 
= {cosep\HA)\VB)+e'^siYiep\VA)\HB)) 
®{cos e m\L a)\L b) + sine m\Ra)\Rb)) ■ 

By varying the crystal phase-matching temperature to change the solid angle of the SPDC 
emission cone, one can adjust the flux ratio between the collinear output mode and the 
non-coUinear output mode, thereby controlling 6m- We set the pump relative phase such 
that the output phase = 0, and the half-wave plate placed in Alice's path in Fig. [T] was 
used to rotate Aice's output polarization by 90° to yield a partially hyperentangled output 
state given by Eq. ([5]). As an input state for Schmidt projection, we adjusted the pump 
relative amplitudes to set Op = 6m = 0- 

B. Implementation of Schmidt projection 
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FIG. 2: (Color online) Entanglement distillation scheme for hyperentangled photons. Alice and 
Bob have the same setup, (a) Schmidt projection transmits only the two relevant terms, \VR) 
and \HL), of the initial state given in Eq. p-cnot gate combines two paths {R, L) into a 
common path for final state extraction and for polarization state analysis, (b) p-cnot gate with 
phase compensation. PBS: polarization beam splitter; HWP: half-wave plate; QWP: quarter-wave 
plate; DP: dove prism. 

We applied single-qubit and two-qubit gates of single-photon two-qubit quantum logic 
to implement Schmidt projection on the partially hyperentangled state of Eq. ([5]), or its 
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expanded form of Eq. ([6]). The SP procedure is to extract from the initial state the maximally 
entangled term \ip)sp = \VR)a\VR)b + \HL)a\HL)b, which is orthogonal to the other two 
terms in Eq. Therefore, \il')sp can be simply isolated by using a polarization beam 
splitter to separate the polarizations and choosing the appropriate beam paths. However, 
the extracted state \i')sp still contains four modes even though V (H) is always associated 
with R (L). It is more useful to convert the extracted state to a more familiar form such as 
the polarization triplet state 



\(j)+) = \HAHB + VAVB)/V2. (9) 

The transformation from \ip)sp to can be accomplished with a polarization controlled 
NOT (p-cnot) gate 2l(] that converts ip)sp to \L), thus bringing the momentum 
qubit to a fixed value, which we can ignore in the final output state \4>~^)- The polarization- 
entangled triplet final state is much easier to analyze using polarization correlation mea- 
surements than dealing with a state containing both L and R paths. 

Fig. [2](a) shows the apparatus for implementing the Schmidt projection method. The 
key step is to project the initial state given by Eq. ([6]) into a subspace composed of \VR) 
and \HL) for both Alice and Bob without destroying the qubits. In principle, this can be 
accomplished by using one PBS in the L path to transmit the if-polarized photon in the 
L path and another PBS in the R path but is rotated by 90° about the propagating axis 
to transmit the ^-polarized photon in that path. The first two terms in Eq. ([6]) are then 
refiected. In the actual setup of Fig. [2] we modify this arrangement by using the HWPi to flip 
the polarization in the R path before PBSi that transmits the desired subspace components 
in both paths. Instead of recovering the polarization state in the R path after PBSi, we 
flip the polarization state in the L path with IIWP2, so that the output after Schmidt 
projection is cos 9 sin e{\HR)A\HR)B + \VL)a\VL)b). Note that HWPi and HWP2 are 



221. We should 



essentially momentum-controlled not (m-cnot) gates described in Ref. 
point out that the arrangement of having one HWP in each path eliminates the need for 
path-length compensation between the two paths that is necessary if only one m-cnot gate 
(i.e., a HWP) is used. 

To eliminate the path dependence of the Schmidt-projected state \HR)a\HR)b + 
\VL)a\VL)b, we fold the L-path output and the i?-path output onto a common path by 
use of a P-CNOT gate, as indicated in Fig. M^a). The flnal output state becomes and 
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we simplify the notation by removing the path designation. The implementation of this P- 
CNOT gate is schematically shown in Fig. [2t^b). The p-cnot gate consisted of a dove prism 
embedded in a polarization Sagnac interferometer that would rotate the incoming image (L 
and R beams) by ±90° depending on the beam polarization 2l|]. The P-CNOT gate has 
the following mapping: \HL) \HU), \HR) \HD), \VL) \VD), \VR) \VU), where 
U (D) refers to the upside (downside) path. For the purpose of relating to the usual |0) 
and |1) states, we identify H, L, and U with |0), and V, R, and D with Our p-cnot 
implementation introduced a fixed amount of phase shifts depending on the input polariza- 
tion ^] and they were compensated by adding a HWP sandwiched by two quarter-wave 
plates (QWPs) whose optic axes were tilted by 45°, as shown in Fig. [2](b). This phase 
compensator was also utilized to add the necessary tt /2 phase shift in our quantum state 
tomography measurements in the next Section. 

To verify that the final output is the maximally-entangled polarization triplet \4>^), we 
analyzed the polarization correlation between Alice's and Bob's coincidence counts. The 
polarization measurements were carried out with polarization analyzers each consisting of 
HWP3, PBS2, and a Si single-photon counter, as shown in Fig.[2]^a). The measurement basis 
for the polarization analyzer was set by HWP3. 



IV. EXPERIMENTAL RESULTS 

Our implementation of the Schmidt projection protocol is best characterized by first 
demonstrating the required input hyp erent angled state of Eq. ([6]) was generated and then 
measuring the degree of entanglement of the output state. Given a general hyp erent angled 
state with arbitrary degrees of entanglement for the polarization and momentum qubit pairs 
of Eq. (EI), 6'p = 9m = 9 is required for the input. In our verification of the input state, we 
characterized the partially hyp erent angled state by comparing the different amounts of en- 
tanglement in the polarization and momentum modes at different 9 = 44°, 41.9°, 39.3°, 35.9°. 
For each 9 the SP protocol was carried out on the input state and we analyzed the extracted 
output state. In particular, we measured how close the extracted state was to the expected 
polarization triplet 
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FIG. 3: Characterization of initial hyperentangled state. Detected coincidences of (a) \L)a\L)b 
and \R)a\R)b terms for momentum input state 9m = 35.9° {6p = 45°), and (b) |1^)a|^)_b and 
\H)a\H)b terms for polarization input state Op = 35.9° {9m = 45°). (c) Initial state of polarization 
qubits tan^^p as a function of initial state of momentum qubits tan^^M- Open squares are the 
measurement results, and the straight line indicates the ideal case. Pump power was 5 mW, and 
detection bandwidth was 1 nm. 

A. Input state characterization 

To characterize the input hyperentangled state of Eq. ([5]), we utilized the SP and p-cnot 
apparatus of Fig. [2t^a) and our ability to control 6m and 6p independently, as described in 
Sect. III. The combination of the SP and p-cnot gate transforms the input state of Eq. ([5]) 
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into 

{cosepsmeM\H)A\H)B + smepcoseM\V)A\V)B)\D)A\D)B . (10) 

To evaluate the distribution of the momentum entanglement, we set 6p = 45° in Eq. ([T0|) 
so that cosOp = sinOp = 1/V2. As discussed in Sect. Ill, we adjusted the relative H and 
V components of the pump to control 9p. We then obtained the ratio of sin^ 9m/ cos^ 9m 
by measuring the coincidence ratio of C{H,H) and C(y,V), where C{X,Y) denotes the 
detected coincidence rate between X-qubit measurement by Alice and F-qubit measurement 
by Bob. Any desired momentum entanglement distribution can be obtained by adjusting 
the SPDC crystal temperature while monitoring this coincidence ratio. Fig. [3]^a) shows the 
measured distribution of the momentum qubits for 9m = 35.9°. For all the coincidence 
measurements in this work, each consisted of an average of 60 1-s measurements without 
background subtraction. 

Similarly, to obtain the distribution of the polarization entanglement, we set 9m = 45° 
by choosing an appropriate PPKTP temperature. After Schmidt projection and the P- 
CNOT gate, we measured the polarization distribution by monitoring the coincidence ratio 
of C{V, V) and C{H, H) to yield tan^ 9 p. Fig.[3]^b) shows the distribution of the polarization 
qubits for the case of 9p = 35.9°. Note that Fig. and (b) have different scales because 
the two measurements were taken at two different temperatures that resulted in different 
output fluxes. More importantly is that Fig. [3]^a) and (b) show that both ratios are nearly 
identical, as required by the Schmidt projection protocol. 

The above measurements for determining the crystal temperature to obtain 9m = 35.9° 
and the pump's H-V component ratio to set 9p = 35.9° were repeated for other 9m and 9p 
values of 44°, 41.9°, and 39.3°. Fig. [21(c) shows the measured polarization coincidence ratio 
tan^ 9p at various measured momentum coincidence ratio tan^ 9m, showing clearly that the 
degree of partial entanglement in the polarization and momentum qubit pairs were well 
matched, 9m ~ 9 p. The ideal distribution ratio is also shown in Fig. [3]^c) as a straight 
line. To obtain the partially hyperentangled state required for the SP protocol (and for the 
measurements in Fig. El^c)), we set 9m = 9p = 9 in the following way. First, we set 6'p = 45° 
that could be confirmed by monitoring the coincidence ratio in the polarization qubits of 
the coUinear i?-path outputs. Then we adjusted the crystal temperature to the value that 
would yield the desired 9m, which we verified with the measurements presented in Fig. [3](a). 
Finally, we set to the same value as 9m = 9 that we confirmed by again measuring the 
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FIG. 4: Coincidence measurements of output state in the H-V basis, (a) Measured coincidence 
rate of the two parahel and two orthogonal polarization terms of the output state for input state 
with 6 = 35.9°. (b) Plot of C{H, H)/C{V,V) of the output state as a function of initial momentum 
qubit ratio tan^ 6m- Pump power was 5 mW, and detection optical bandwidth was 1 nm. 

polarization coincidence ratio of the collinear i?-path outputs. 



B. Characterization of Schmidt projected output state 

After the distillation procedure, we expect the final state to be a maximally polarization- 
entangled triplet state of Eq. We first analyzed the output state by measuring 
the coincidences between Alice and Bob in the H-V basis, that for should yield equal 
probability of detecting parallel polarizations but zero probability of detecting orthogonal 
polarizations. Fig. Hl^a) shows the four measured coincidence rates HaHb, H^Vb, VaHa, 
and VaVb for an input state with 6 ~ 35.9°. The measurement results clearly show that the 
output state had nearly balanced \H)a\H)b and terms and negligible \H)a\V)b 

and \V)a\H)b components. We then made the same coincidence measurements for input 
states with different values of 6. The coincidence rates for orthogonal polarizations HaVb 
and VaHb were negligible, as expected for \4>'^). Fig. |l](b) plots the ratio of coincidence rates 
C{H, H)/C(y, V) of the output state (solid triangles) as a function of the input state set by 
6. The measured coincidence ratios lie close to the expected value of unity for a maximally 
entangled output state with less than 5% error. 

To show that the distilled output was indeed entangled as a coherent superposition, and 
not as a classical mixture, of \H)a\H)b and we made additional measurements 



for verification. We first performed quantum state tomography |3l| on the output state. 
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FIG. 5: (a) Real part, and (b) imaginary part of measured density matrix of the Schmidt projected 
output state for initial state with 6 = 35.9°. 
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FIG. 6: Measured quantum-interference visibilities in the ±45° diagonal basis {\H) it \ V)). 

We can reconstruct the density matrix of Alice's and Bob's polarization qubits by making 
coincidence measurements for the 16 combinations of their polarization analyzers set along 
\H),\V), {\H) + \V))/V2, and (|i/)-i|1/))/v/2. The real and imaginary parts of the output 
density matrix p for the case 9 = 35.9° are displayed in Fig. [5](a) and (b), showing clearly 
that the state is with a fidelity, (0+|p|0+), of 0.952. 

Another useful indicator of the output state coherence is two-photon quantum interference 
visibility measured in the ±45° diagonal basis {\H) ±\V)). Quantum interference visibility is 
defined as {Cmax — Cmin)/{Cmax + Cmin), whcrc Cmax IS the maximum coincidence counts and 



15 




0.0 0.5 1.0 1.5 2.0 

Entanglement in initial state (E(i/r)) 



FIG. 7: Measured efficiency of Schmidt projection implementation. Expected number of maximally 
entangled pairs after Schmidt projection is applied to a pair of hyperentangled state as a function 
of entanglement in the initial state {E{%1))). Open circle is the measured number and the solid 
curve is the theoretical prediction when Schmidt projection is applied to a hyperentangled state. 
The dashed curve shows the maximum number of maximally entangle pairs obtainable by Schmidt 
projection when the input consists of an infinite number of pairs. 

Cmin is the minimum coincidence counts. Fig. [6] plots the measured visibilities for different 
input states, all showing ~90% visibility except one lower-visibility point that was caused by 
a slight misalignment of the apparatus. The visibilities in Fig. [6] were primarily limited by 
the imperfect p-cnot gate that had a classical visibility of only ~93% [21']. Another possible 
source of the visibility loss is the slight spectral mismatch of the output photons from the L 
and R output spatial modes. Even though the spectra of the output photons were mainly 
determined by the two interference filters, a slight mismatch in the filtered output spectra 
could happen because slightly different emission angles for the L and R paths (see Fig. [1]) 
would yield slightly different spectral distribution. 

From the different state analysis measurement results shown in Fig. HI [5l and El we 
conclude that our distillation outputs were in the state independent of the input states 
with different 9 values. One of the key advantages of the SP protocol is its high distillation 
efficiency, especially for a large number of initial pairs. We have measured the efficiency 
of our Schmidt projection implementation in terms of the expected number of maximally 
entangled pairs remaining in the final state as we varied the amount of initial entanglement 
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E{iIj) that we plot in Fig. [7] as open circles. We find that the measured efficiency is in 
good agreement with the theoretically calculated values (solid curve). For comparison, the 
dashed line in Fig. [7] shows the maximum number of maximally entangled pairs obtainable 
by the Schmidt projection protocol when it is applied to an infinite number of in put pairs. 
In general, Schmidt projection shows a lower yield than Procrustean methods [zl. I12I when 
it is applied to less than 5 pairs, as in our case. However, even with a small number of initial 
pairs, Schmidt projection still has the practical advantage that it is not necessary to adjust 
the distillation setup as a function of the input states. 



V. DISCUSSION 



In this work, we have experimentally demonstrated the Schmidt projection method that 
is considered one of the most powerful distillation protocols. We employed hyperentan- 
gled photons produced by spontaneous parametric down-conversion to generate pairs of 
partially-entangled states. We developed a technique to independently control the amount 
of entanglement in the polarization (by pump adjustments) and momentum (by crystal tem- 
perature) degrees of freedom of the hyperentangled photon pairs. This entanglement control 
capability was utilized to prepare the identical pairs of partially entangled qubits as the ini- 
tial state for the SP protocol. Schmidt projection on the hyperentangled photon pairs was 
realized by a PBS and two m-cnot gates, and we used a p-cnot gate to convert the output 
into a polarization Bell state. Our experimental results show that Schmidt projection can 
distill a pair of maximally entangled qubits in the form of a polarization-entangled triplet 
and it is independent of the degree of entanglement in the initial state. Furthermore, the 
measured distillation efficiencies agree with the theoretical prediction. We note that Schmidt 
projection can be selectively applied to perform entanglement purification. If classical com- 
munication is used to compare the projection results, then certain types of mixed states can 
be purified, such as p = T,p\i/j) {ipl + (1 — p)|01)(01|, as well as a mixture of two identically 
decohered pairs [14 1. 

It is important to recognize that the key to our Schmidt projection implementation is that 
distinct tensor products of local qubits, |0)i|0)2, |0)i|l)2, |l)i|0)2, and |l)i|l)2, are encoded 
in physically distinguishable states, \VL), \VR), \HL), \HR). Therefore the same concept 
can be applied to other qubit systems such as cavity quantum electrodynamics (cQED) 
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FIG. 8: Schematic illustration of applying Schmidt projection to cQED. (a) Alice and Bob have 
two pairs of entangled atoms (Ai-Bi, A2-B2) trapped in an optical cavity. Pulses from two lasers, 



Li and L2, are applied in a specific sequence to implement adiabatic transfer of the dark state 



(b) Proposed electronic-level structure for implementing the SP protocol. Atoms Ai and A2 are of 



the same kind (e.g., Rb [3^), and different levels are used to encode the initial qubit state. After 
adiabatic passage, Ai will be in the |0)i state, and A2 will be a superposition of all four states 
(10)2-13)2). 



32, l33|] or trapped ion [3J] systems, where multiple-qubit product states can be mapped into 
different internal states within the same atom. For example, in a cQED system illustrated 
in Fig. [Hl^a), Alice has two trapped atoms each of which is entangled with the corresponding 
trapped atom in Bob's cavity, where the qubits are stored in the ground-state Zeeman levels. 
Consider the case in which the first qubit is stored in |0)i and |l)i of Ai (Bi) and the second 
qubit is in |0)2 and |1)2 of A2 (B2) on Alice's (Bob's) side as shown in Fig. [8](a). One can 



coherently transfer the quantum states stored in Alice's two atoms completely to A2 whi 



e 



leaving Ai in |0)i by adiabatic passage via a dark state of the two-atom + cavity system 32|. 
That is, \p)i\q)2 => |0)i|p5')2 when is in a binary representation. By mapping all four 
possible quantum states into four internal states of A2, the undesirable states such as |0)2 
or 13)2 can be eliminated by driving the cycling transitions from the undesirable states to 
an auxiliary state (not shown in the level structure) 3^, and the process can be confirmed 
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by the induced fluorescence. If no fluorescence is observed, we can conclude that Schmidt 
projection has successfully taken place. 

In the trapped- ion case, one can use the quantized center-of-mass vibrational mode ('bus 
mode') in a similar way as in the cQED case. When there are two ions, the electronic state 



of the flrst ion can be transferred to the bus mode [3^ with a red-sideband transition, and 
this bus mode can be transferred to the second ion while mapping four different two-qubit 
states to four distinct internal states of the second ion. The rest of the procedure is similar 
to the cQED case. We believe that further studies of Schmidt-projection implementations 
in atoms and ions and for mixed states should enhance our expanding collection of tools in 
quantum information science. 

This work was partially supported by the Hewlett Packard-MIT Alliance. 
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